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ABSTRACT 



The Cepheus E outflow has been studied in the mid and far infrared using the 
ISO CAM and LWS instruments, and at miUimetric wavelengths using OVRO. In 
the near and mid-IR, its morphology is similar to that expected for a jet driven 
outflow, where the leading bow shocks entrain and accelerate the surrounding 
molecular gas. As expected, fine structure atomic/ionic emission lines arise from 
the bow shocks, at both the Mach Disk and the stagnation tip, where J-shocks 
are dominant. The H2, H2O and CO molecular emission could arise further 
'downstream' at the bow shock wings where the shocks (v = 8 — 35 km s~^) are 
oblique and more hkely to be C-type. The ^^CO emission arises from entrained 
molecular gas and a compact high velocity emission is observed, together with 
an extended low velocity component that almost coincides spatially with the H2 
near-IR emission. The miUimetric continuum emission shows two sources. We 
identify one of them with IRAS 23011+6126, postulating is the driver of the 
Cepheus E outflow; the other, also an embedded source, is hkely to be driving one 
of other outflows observed in the region. Finally, we suggest that the strong [C II] 
158 /im emission must originate from an extended photo-dissociation region, very 
hkely excited by the nearby Cepheus 0B3 association. 



Subject headings: stars: formation — stars: pre-main-sequence — ISM: jets and 
outflows — infrared: ISM: lines and bands 
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1. Introduction 

Low mass proto-stars produce outflows that interact mostly through shocks with the 
surrounding medium (e. g. Reipurth & Raga |1999| ). Nevertheless, the details of the earliest 



evolutionary steps remain somewhat unclear. The embedded flow at Cepheus E (hereafter 
Cep E) displays all the characteristics of an extremely young object still surrounded by a thick 
cocoon of cold gas (e. g. Cernicharo et al. pOOOD , and hence provides an ideal opportunity to 



study the development of molecular /ionic outflows and their surrounding environments during 
the earliest stages of star formation. 



Cep E was identified originally as a cloud core in molecular CO observations (Sargent 1977|) 



and subsequently, an associated outflow was discovered (see Fukui |1989| ; Hodapp |1994| ). The 



region around Cep E has proved more complicated than the early CO and the K' observations 
indicated; there are at least two other outflows besides the original compact 'North-South' 
H2 flow identified by Ladd & Hodapp (|1997D . One is nearly perpendicular to the main H2 
flow and appears as a faint chain of H2 1-0 S(l) 2.121 fim emission knots towards the West 
(Eisloffel et al. |1996D ; the other is seen in the CO J = 1-0 line, is extended by ~ 90" along 
a PA ~ -45° (Ladd & Howe 1997| ). All three flows arise within a 5" region around IRAS 



23011+6126, which presumably encompasses the driving source (Lefloch, Eisloffel & Lazareff 
1996| ; Noriega- Crespo, Garnavich & Molinari |1998| , hereafter NCGM; Cernicharo et al. 



pOOOp . For a low mass object IRAS 23011+6126 is estimated to have a very high bolometric 



luminosity, ~ 70 L© (Lefloch et al. |1996|) , suggesting that we are dealing with a relatively 
massive protostar, or with a low-mass protostar in an unusually early evolutionary phase and 
considerable higher accretion luminosity. 

The blueshifted South lobe of the original H2 Cep E outflow is observed at optical 
wavelengths as Herbig-Haro HH 337, in the shock excited Ha and [SII] AA6717/31 emission 
lines (Noriega- Crespo |1997| ; Devine, Reipurth & Bally |1997| ). This collisionally excited 
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emission is characteristic of an extremely low excitation HH object (~ 1% ionization) with an 
anomalously high electron density, ~ 4100 cm~^, leading to an unusually high pre-shock 
density of ~ 10^ cm~'^ (Ayala et al. pOOOj ; hereafter AYA). Both the Ha and [SII] emission 



can be modeled using a J shock of 'm shock ~ 15-30 km s ^ (AYA), although the near-IR H2 
spectra agree better with C-shock models with \shock ~ 35 km s^^ (Ladd & Hodapp |1997|) 



It has been suggested that a Mach disk produces the optical emission while the IR emission 
arises from a bow shock (AYA). 

Some of the first three dimensional simulations of molecular jets were intended to model 



the complex structure of the H2 Cep E outflow (Suttner et al. |1997| ). Surprisingly they 
demanded relatively large (> lO^cm"^) densities, which could be a sign of a very early stage 
of jet development. This interpretation is consistent with the estimated dynamical age of the 
Cep E outflow, ~ 3 X 10^ years (Eisloffel et al. 1996 ), assuming a distance of 730 pc. The 
above argument sets a scene in which Cep E is a very young outflow breaking through its 
placental molecular core. 

Our new observations from the Infrared Space Observatory (ISOQ) and the Owens 
Valley Radio Observatory (OVRO) millimeter wave array at far/mid infrared and millimetric 
wavelengths at higher sensitivity and spatial resolution, permit not only a probe of deeper 
regions of the molecular environment surrounding the outflow, but also a determination of 
the main excitation and cooling mechanisms. Both are essential to understand the dynamical 
evolution of very young proto-stellar objects. 

In § 2 and § 3 we describe the ISO and OVRO observations and data reduction procedures, 
respectively. The infrared results are presented in § 4 and further analyzed in § 5, where we 



^ISO is an ESA project with the instruments funded by ESA Members States (especially 
the PI countries: France, Germany, the Netherlands and the United Kingdom) and with the 
participation of ISAS and NASA. 
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derive the gas physical parameters (such as temperature, density, column density and size 
of the emitting region) and study the origin of this emission through J and C-shocks and 
photo-dissociation regions (PDRs). In § 6 we analyze the millimeter results, using the ^^CO 
line emission to derive such physical parameters of the outflow as dynamical time scale, mass 
and mass loss rate, and continuum observations to provide an alternate estimate of the mass 
of the source. The source spectral energy distribution (SED) is examined in § 7 and in § 8 we 
summarize our results. 



ISO Observations and data reduction 



The Cep E outflow has been observed with two ISO instruments, the Long Wavelength 



Spectrometer (LWS, Clegg et al. |1996| ) and the infrared camera (ISOCAM, Cesarsky et 
al. [L996| ). The observations were taken during revolutions 566 (LWS: TDT 56601113 and 
56600912) and 792 (CAM: TDT 79200740). The LWSOl grating mode was used to acquire 
low resolution (R~ 200) spectra from 43 to 197 /im at two different locations along the 
flow (a(2000) = 23h03ml3.11s, 5(2000) = 6r42'59.5" for the North lobe and a(2000) = 
23h03ml2.73s, 5(2000) = 61°41'56.5" for the South lobe). Figure 1 shows the two pointing 
positions. These spectra are made up of 11 full grating scans oversampled at 1/4 of a 
resolution element (equivalent to ~ 0.07 fim for A<90 ^m and to ~ 0.15 fim for A>90 fim ). 
The integration time per spectral element was 22 s, corresponding to a total integration time 
of 1890 s per pointing. The LWS beam is assumed to be ~ 80" in diameter (ISO Handbook, 



IV, 4.3.2). The ISOCAM (Cesarsky et al. |1996| ) observations were taken using the Circular 



Variable Filter (CVF) with a total FOV of 3' x3', and pixel FOV of 6" centered on IRAS 
23011+6126 (a;(2000) = 23h03ml3.02s, 5(2000) = 6r42'26.5"). The CVF spectra covered the 
5.01 — 16.77 yum wavelength range with a resolving power of ~ 40. Each wavelength was 
observed for ~ 23 sec, for a total on source observing time of 3550 sec. 
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The LWS Interactive Analysis [| (LIA) Package, Version 7.2 was used to correct the 
detector dark currents and drifts. Corrected and scan averaged spectra were then analyzed 
using the ISO Spectral Analysis Package^ (ISAP) Version 2.0. Likewise, the CAM-CVF data 
were reduced using the Cam Interactive Analysis (CIA) standard procedures, and spectra 
from the individual pixels were extracted using XCVF (Ali pOOO|) , and further analyzed with 
ISAP. 

3. Millimeter observations and data reduction 

High resolution millimeter interferometric observations were obtained during 1997 and 
2000 using the OVRO millimeter-wave array. The array comprises six 10.4 m antennas that 
were deployed in three different configurations offering baselines extending from the shadowing 
limit to 220 m. To ensure the best possible sensitivity to the extended outfiow emission, the 
data from the most compact available configuration is included. Each telescope is equipped 
with cryogenically cooled SIS receivers providing typical system temperatures of ~ 300 K at 
110 GHz and ~ 1200 K at 222 GHz. Continuum observations employed an analog correlator 
for total bandwidths of ~ 2 GHz and 1 GHz at 222 GHz and 110 GHz (1.3 and 2.7 mm) 
respectively. The digital correlator was configured to observe simultaneously the ^^CO(l-O) 
and C^^O(l-O) transitions at 110.2 and 109.8 GHz respectively. For the former isotopomer, 
spectral resolution was 1.3 km s^^ over a ~ 80 km s^^ bandwidth, and for the latter, it was 
0.3 km s~^ over ~ 20 km s^^. The half-power beamwidth is ~ 65" at 110 GHz, and ~ 33" at 
222 GHz. 



^LIA is available at |http: / / www.ipac.caltech.edu / iso / Iws/lia / lia.html 
^ISAP is available at Pittp:/ /www. ipac.caltech.edu/iso/isap/isap. html 
^CIA is available at Pittp: / /www. ipac.caltech.edu/iso/ cam/ camtools.html 
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Complex gain calibration was ensured by frequent observations of the quasar BG 0059+58 
(= [MFN99] 0059+581). The typical observing cycle alternated scans on Cep E and BG 
0059+58 with a period of ~ 20-30 min. 3C 345 and/or 3C 84 were used as passband calibrators 
while scans on Uranus and/or Neptune were used to set the flux scale. The expected flux 
calibration uncertainty is ~ 15% at 2.7 mm, and 20-30% at 1.3 mm. Data calibration and 
editing have been performed using the OVRO-MMA software package (Scoville et al. [19931 ). 
Imaging and analysis of the calibrated {u, v) datasets were carried out with the NRAO-AIPS 
package. All the maps presented in this paper have been obtained with the AIPS IMAGR 
task using natural weighting of the {u, v) data for the molecular outflow maps and the low 
spatial resolution continuum maps, while "ROBUST=0" weighting was used for the 1.3 mm 
high resolution map. The beam sizes and noise levels for the OVRO observations are shown in 
Table 1. 

4. ISO Results 
4.1. ISO LWS 

Figure 2 shows the LWS far-infrared spectra obtained at the two pointing positions. 
The strongest features are the atomic fine structure lines of [O I] 63 /im and [G II] 158 fim. 
However, the spectra are also rich in molecular lines from the rotational transitions of GO and 
H2O (in both ortho and para form). Two OH lines are also detected in the spectra, but they 
are too faint to be included in the analysis. 

Integrated line fluxes, measured using a gaussian fit, are shown in Table 2, with 
1-0" uncertainties given in parenthesis. In some cases the FWHM was kept fixed to the 
corresponding instrumental resolution element (i.e. 0.29 ^m and 0.6 fim for the SW and LW 
detectors, respectively). This technique was particularly important in the case of faint or /and 
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blended lines. The uncertainties were obtained either by using ISAP estimated values from 
the gaussian fit, or by adopting the RMS of the fitted baseline times the resolution element. 
As expected, both values are very similar. 



4.2. ISOCAM CVF 

CAM CVF provides spectra covering the 5-17 fim range over a 32x32 pixel area. Each 
plane of the cube corresponds to a different A. Figure 1 shows, as was previously noticed by 
NCMG, how the morphology at 2 and 7 fim in the H2 emission is quite similar. A more 
detailed description of the CVF emission at each pixel is given in Figure 3. Except at the 
source position (pixel (16,16)), the 5 to 16 fim spectra are dominated by the H2 ground 
rotational lines from S(7) to S(2). Integrated line fluxes can be found in Table 3. These are 
corrected for extinction, using an estimated E(B-V) = 1.0 (Lefioch et al. 1996| ; AYA) and the 



extinction curve from Draine (|1989| ). The central source is clearly detected and displays strong 
silicate absorption at 9-12 fim and CO2 ice at 15.2 fim and a number of ice bands (see also 
Cernicharo et al. |2000|) . Figure 4 displays the images corresponding to the H2 0-0 S(2) to S(7) 



lines. In accord with expectations based on ground based optical and H2 observations at 2 /im, 
the south lobe is brighter at shorter wavelengths, as the flow breaks through the cloud. The 
north lobe, buried deeper in the cloud, has fainter S(3) 9.66 //m and S(2) 12.28 fim emission, 
due probably to the broad second maxima at 9.7 /xm due to silicates (Si-0) stretch mode in 
the extinction curve (Draine |1989|, Mathis |1990|). 



5. Discussion 



5.1. The Gas Physical Parameters 
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5.1.1. CO and H2O Emission 

The CO and H2O emission can be studied by means of a Large Velocity Gradient (LVG) 
model in spherical geometry, taking into account the first 35 rotational levels of CO and the 
first 30 rotational levels of 0-H2O. The collisional rates of Schinke et al. (|1985|) were used for 
CO and those of Green, Maluendes and Mc Lean ( |1993 ) and Phillips et al. ( |199(j| ) for H2O. 



There are six main parameters to fit: temperature (T), H2 density (n(H2)), CO and 0-H2O 
column densities (N(CO) and N(o-H20)), intrinsic width of the line (ve^p) and angular size of 
the emitting area {Q). To determine all of them simultaneously, two initial assumptions are 
necessary: (1) the CO and 0-H2O emission arise from the same region; (2) the CO lines are 
optically thin, i.e. their line fluxes do not depend on their intrinsic widths, while the 0-H2O 
lines are optically thick. Since the LWS resolution is too poor to establish size and position of 
the emitting region, it is difficult to justify (1) from our data. However, in a 15" beam, the 
H2O emission in HH 7-11 appears to arise from the shocked regions already seen in the CO 
emission bullets (Cernicharo, Gonzalez- Alfonso & Bachiller ( |1996[ ) and Bachiller & Cernicharo 



1990| ). The adoption of assumption (2) will be confirmed by our models. Figure 5 shows the 



CO line fiuxes as a function of the rotational quantum number, J„p, for Cep E North and 
South lobes. The shape of the distribution is determined only by T and n(H2), while the 
absolute fiuxes are determined by T, n(H2) and fixN(CO), which is proportional to the total 
number of CO molecules within the emitting volume. 

Two extreme cases for each lobe are consistent with the observations. All the models 
predict that the CO lines are effectively optically thin. Figure 5 shows that the agreement is 
good for Jup<22; but for higher J up the models underestimate the observed fiuxes. This could 
indicate the existence of a hotter CO gas component as suggested by Spinoglio et al. (|2000| ) 
for T Tau. 



Assuming that the CO and the 0-H2O line emission is from the same gas, the temperatures 
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and densities derived from the CO lines can be used to model the 0-H2O line fluxes. Since 
these lines are optically thick, their absolute fluxes depend on all the free parameters in the 
model, T, n(H2), N(o-H20), Vg^p and Q. Their relative fluxes (or line ratios), on the other 
hand, do not depend on Q. Thus, with T and n(H2) known, the line ratios were modeled and 
the values of N(o-H20) and Vexp were constrained. Models of the absolute fluxes then allowed 
a derivation of Q. The 0-H2O lines at 179, 174, 113 and 75 /im were used because they have 
the highest signal to noise ratio and because the 0-H2O 179 /im/174 yum ratio is very sensitive 
to T and n(H2). 

A summary of the LVG model parameters consistent with the observations can be found 
in Table 4. A comparison between the LVG models for 0-H2O and the observations is shown 
in Figure 5. The diameter range estimated for the 0-H2O emitting region is comparable to 
the extension of the vibrational H2 1-0 S(l) 2.12 ^m emission, 2700-6900 AU (3.7" - 9.5") at 
a distance of 730 pc. Knowing the emitting area, we calculate N(CO) and hence the H2O/H2 
abundance ratio, assuming a standard value of C0/H2~ 10"'^ in shock regions. For H2O, 
an ortho/para ratio of 3 was adopted, leading to predicted P-H2O fluxes within 3-a of the 
observed values. The low temperature models give H20/H2~ 4.0x10"^ and 1.5x10^^ for the 
North and South lobe respectively. The high temperature models give H20/H2~ 2.5x10"^ 
(North) and 2.4x10^^ (South). Taking into account these estimates together with their 
uncertainties, it can be concluded that H2O abundances are between 100 and 4000 times 
greater than expected for the ambient medium in cold quiescent conditions, x(H20)~ 10^'' 
(Bergin et al. |1998| ; Moneti, Cernicharo & Pardo |2000|) , a result commonly found in young 



stellar objects (Liseau et al. p.996| ; Cernicharo et al. |1996| ; Harwit et al. |1998| ; Gonzalez-Alfonso 



et al. |1998| ; Ceccarelli et al. |1999| ; Spinoglio et al. |2000D . This significant enhancement in the 



water abundance within outflow regions can be attributed to the vaporization of volatile grain 
mantles and the increased rate coefficients for specific gas-phase reactions at the elevated gas 
temperatures reached behind the shocks, and it is expected to persist for ~ 10^ years (Bergin 
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et al. [1998D . 



According to the LVG models, the high-J CO emission detected by ISO-LWS indicates a 
temperature between 220 K and 1200 K. Since the contribution to low-J emission from gas at 
such temperatures is neghgible, the ISO observations must be tracing only the hottest shocked 
molecular material. By contrast, millimetric observations of the low-J emission (J=2-l, Lefloch 
et al. |1996| and J=4-3, Hatchell et al. |1999|) provide information about the entrained gas. 



5.1.2. H2 Emission 

The pure rotational H2 lines seen in the ISOCAM CVF spectra arise from quadrupole 
transitions that are always optically thin. Due to their low Einstein coefficients, thermalization 
is readily possible for moderate volume densities. Hence, LTE calculations are precise enough 
to compute H2 column densities. In outflows from YSOs, the H2 emission is mostly produced 
by coUisional excitation driven by shocks, but to gain some insight into the overall nature of 
the gas, local thermodynamic equilibrium is commonly assumed. It is then possible to estimate 
the excitation temperature (Tex) and column density of the H2 gas using ln(Nj/gj) = ln(Ntoi 
/Q) - Ej/(k Tex), where = 4 tt Fj/{Q Ej gj Aj). As usual, Aj is the transition probability 
(Turner, Kirby-Docken & Dalgarno [19771) and E^ the energy of the upper level (in cgs units) 
(Dabrowski |1984| ). The dereddened fluxes, Fj/Q, are listed in Table 3. The partition function 
(Q) was calculated assuming an ortho-para ratio of 3. The excitation diagrams, ln(Nj/gj) vs 

and the derived excitation temperatures and column densities for several pixels along the 
flow are shown in Figure 6. Given Nj and Tex, a model spectrum can be created for all H2 0-0 
lines. A comparison one of these models with the observations is shown in Figure 6. A single 
component model predicts the intensities from the S(3) to the S(7) lines very well, but the 
intensity of S(2) is systematically underestimated. 



- 13 - 



5.2. Interpreting the Shock Excited Emission 

5.2.1. C-shocks 

The ground rotational H2 emission observed in YSOs is mostly attributed to collisional 
excitation from shocks, either C-type or J-type (Draine & McKee |1993| ). In a nutshell, 



C-type shocks are inherently magnetic and less harsh on molecules because Alfvenic waves 
prepare the gas to the hydrodynamic shock. J-type shocks are mostly hydrodynamic and can 
dissociate H2 molecules at lower shock velocities. In the mid and far infrared, where molecular 
transitions are detected, C-type shocks are expected to be more important. This is confirmed 
by the comparison of the CO and II2 pure rotational line fluxes predicted by C-shocks models 
for Vshock ~ 10 — 40 km s~^, and gas densities 10"^ — 10^ cm~^ (Kaufman & Neufeld |1996|) , and 
J-shocks models for y shock ~ 30 — 150 km s~^ and the same gas density range (Hollenbach & 



McKee |1989| ). C-shock fluxes are at least 10 times brighter than the J-shocks (and this factor 
can be much larger depending on the shock conditions). A similar situation is expected for 
H2O, and the following analysis assumes that only C-shocks contribute. 

The shock velocity and pre-shock density are estimated using the pure rotational H2 line 
fluxes and their ratios. From a comparison of these values with C-shock models (Kaufman & 
Neufeld |1996| ), we determine a range of 20—40 km s~^ shock velocities, but over a wide range 



of pre-shock gas densities. The Vghock value of ~ 35 km s ^ derived from the near infrared H2 



emission lines (Ladd & Hodapp 1997 ), falls within this range. An estimate of the pre-shock 
gas density can be obtained with the additional assumption that the emission from the pure 
rotational H2 and CO lines arise from the same region. If so, the CO/H2 cooling ratio can 
be used in conjunction with the shock velocities. Since this ratio encompasses all cooling, 
the fraction of lines lying outside the LWS and CVF wavelength coverage must be taken into 
account. The CO, 0-H2O and P-II2O cooling rates derived from the LVG models can be 
found in Table 4. For CO, the LVG models neglect the contribution from lines with J^p >22, 
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which may belong to a hotter component. The H2 coohng was estimated using the simple 
LTE model and, after taking into account the difference between the LWS and CAM beam, 
is 1.6x10^^^ W cm~^ for both the North and South lobes. Figure 7 shows the comparison 
between the corrected CO/H2 cooling ratios and those predicted by C-type shock models 
(Kaufman & Neufeld |1996|) . Within the Vshock range of 20—40 km s~^, the models predict a 
pre-shock density of 10^ — 2.5 x 10^ cm~'^. 

Figure 8 shows the predicted CO/H2O cooling ratio as a function of pre-shock density and 
shock velocity for C-shocks models (Kaufman & Neufeld |1996| ). The observed (and corrected) 



ratio for the Cep E outflow indicates shock velocities of 8-14 km s~^, inconsistent with the 
values derived from the CO/H2 ratios. This could be due to the fact that the brightest 
emission regions in H2 and CO do not seem to coexist spatially, as has been assumed, and 
therefore trace different gas conditions. Indeed, ^^CO presented in Figure 9 suggest this may 
be the case. We discuss this further below. The spatial resolution available with LWS and 
CAM spectra prevent further discrimination between shock diagnostics. 



5.2.2. J-shocks and Photo-Dissociation Regions 

The integrated fluxes of the [O I] 63 //m and [C II] 158 /im emission are comparable for 
both the North and South outflows (see Table 2). Since the strength of [C II] line predicted by 
shock models is several orders of magnitude fainter than the [O I] line, a fraction of the cooling 
is due to the presence of a PDR or photodissociation region (Hollenbach & Tielens p,999| ). Most 



of the [O I] 63 fim emission is expected to be produced by shock excitation, but a percentage 
comes from a PDR. And so, to best use the [O I] 63 /im and [C II] 158 /im lines as diagnostics 
of the PDR parameters, an estimate of the [O I] collisional fraction is necessary. For this, 
we use the additional information provided by the optical forbidden atomic lines (e. g. [O I] 
A6300, [S II] AA 6717/31), and that they are generated by J-shocks, very likely at the Mach 
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Disk (AYA). Implicitly we are assuming that the molecular emission arises at the bow shock 
wings, where C-shocks dominate the emission. 

Near infrared and optical observations suggest a range of J-type shock velocities of 
15 — 35 km s^^. Assuming Yghock = 30 km s^\ the predicted ratio of [O I] 63 /im to [O I] 6300 
A is ~ 10 for n = 10^ cm'^ (HoUenbach & McKee |1989| ). Since the observed [O I] A6300 
flux from the south lobe is 4.7 x 10"^^ erg s~^ cm~^ arcsec"^, the predicted [O I] 63 /im flux 
over the entire object (r ~ 2.5 ") is 9.2 x 10~^^ erg s~^ cm~^. This is of order ~ 25% lower 
than the LWS flux of 12 x 10^^^ erg s~^ cm~^. From a comparison of J-shock and C-shock 
models (c.f. Hollenbach & Mckee |1989| , Draine, Roberge & Dalgarno |1983| ; Timmermann 



1998| ), we estimate that at n = 10^ cm ^ the contribution of C-shocks to the [O I] 63 /im flux 



is about 5-10% that of J-shocks. We will assume that about ~ 20% of the [O I] 63 //m flux 
corresponds to the PDR component. However, the percentage drops as the shock velocity 
decreases, and indeed for y shock = 20 km s~^, the coUisional emission matches the observed 
LWS flux. The corrected ratio of [O I] 63 fim to [C II] 158 /xm emission is then 0.30. 
With the corrected value of [O I] 63 /im/[C II] 158 /xm , and the observed LWS [C II] flux 
of 6.8x10"^ erg s~^ cm~^ sr~^, PDR models can be used to determine the far-ultraviolet 
(FUV) flux and gas densities required to reproduce the observations. Two of the models of 
Kaufman et al.( |1999| ) satisfy the above constrains. For one, the incident FUV flux in units of 
1.6 X 10~'^erg s^^ cm~^. Go, is 30, at a density of 60 cm~^, while for the other Go = 6, at a 
density of 5500 cm~^. 

The integrated flux in the PAHs bands (e.g. at 6.2, 7.7 and 11.3 /im) is correlated with 
the FUV radiation fleld (Boulanger et al. |1998|) . Thus the existence of a global FUV fleld can 
be indirectly inferred from the presence of extended PAHs emission over the entire FOV of the 
CVF ISOCAM observations (see Figure 3). However, the PAHs emission could arise anywhere 
along the line of sight to Cep E. 
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Cep E is near the Cepheus 0B3 association (Sargent 1977 ) with members at distances 
of 500 - 1000 pc (Crawford & Barnes |1970| Jordi, Trullois & Galadi-Enn'quez |1996|) . Since 
there are at least 42 bona fide O and B stars in the association (Garmany |1973| ; Naylor & 



Fabian |1999|) , the diffuse FUV field should be sufficient to generate the observed [C II] ffux. 
One member, HD 217086 (07n), could by itself account for a FUV ffux with Go ~ 6 if it 
is at 725 pc (Garmany |1973| ), similar to the distance to Cep E, and separated from Cep E 
by ~ 80'. Hipparcos measurements of HD 217086 give a parallax of 1.20± 0.92 milli-arcsec, 
approximately 100 pc farther away, but the errors are large. 



6. Results from the Millimetric Observations 

6.1. i^CO Emission 

The Cep E outffow can clearly be seen in the ^^CO J = 1-0 transition at 110.2 GHz, 
shown in Figure 9. One of the most interesting features of this ffgure is how the ^^CO emission 
is almost bounded by the H2 outffow lobes, except at the edge of the south ffow, which is 
breaking through the molecular cloud. Similar behavior is observed for the high velocity 
CO gas in the HH 211 embedded outffow (Gueth & Guilloteau Fig 4), and in HH 1/2 

(Moro-Martin et al. M^, Fig 2c). 



In LTE approximation, using an excitation temperature of 20 K (Ladd & Hodapp 
1997] ), an abundance of ^^C0/H2 ~ 1.2x10^® and an inclination angle of 45°, the following 



parameters are derived for the red and blue-shifted outffow lobes respectively: dynamical time 
scales of 4000 and 8000 years, masses of 0.08 and 0.05 M©, kinetic luminosities of 0.1 L© for 
both, and mass loss rates of 7x10"^ and 6x10^^ M0/year. Nevertheless, these estimates are 
quite uncertain due to spatial ffltering problems, core contamination, and the fact that the 
edges of the outffow are more than 20 " from the phase center of the map, and thus affected 
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by primary beam attenuation. 

Figure 10 shows a position-velocity diagram along the outflow axis (PA=10°) with the 
ambient cloud at -11 km s^^and the outflow extending approximately from to -20 km s~^. 
Faint (2-cr) emission is also seen at ~ 30 and -40 km s~^. We may be tracing high velocity 
gas that is concentrated nearer to the source than its low velocity counterpart. The left panel 
of Figure 11 shows the emission from this gas integrated from -47 to -20 km s^^ and from 
to 25 km s^^. There is no correspondence between this high velocity ^^CO component and 
the ^2 CO 4-3 bullets detected by Hatchell, Fuller & Ladd ( |I999| ). The latter are ~ 10" away 
from the core and have velocities of -120 and 60 km s~^ LSR (a velocity range not covered by 
our observations). A position- velocity diagram along the axis of the second outflow detected 
in ^^CO by Ladd & Hodapp (|1997| ) (PA=-45°) was also constructed but no emission was 
detected. 

Images of the C^^O J = 1-0 emission at 109.782182 GHz were also made, and show a core 
surrounded by a poorly imaged halo. There is no evidence of the Cep E outflow, which is not 
unexpected from this optically thin emission. 

6.2. Continuum Emission at 222 GHz and 110 GHz 

The right panels of Figure 11 display the continuum emission at 222 GHz and 110 
GHz around IRAS 23011+6126. The 222 GHz image shows two unresolved sources in a 
1.38" X 1.06" beam, with integrated fluxes of 90 mJy and 67 mJy respectively. At 110 GHz, 
the continuum emission is resolved and has a diameter of ~ 3" or ~ 2200 AU (at 730 pc) and 
a flux of 35 mJy. It is tempting to interpret this observations as a double source surrounded by 
an envelope. The presence of a second source strongly suggests that the multiple outflows are 
due to multiple sources and not to precession. The projected distance between these sources is 
of ~ 1.4" or ~ 1000 AU at 730 pc. This separation implies a long orbital period, even for an 
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intermediate mass system, with a range of 7 x 10^ - 10"* yrs, i.e. larger than the kinematical 
age of the outflow. 



Using a dust temperature of 18 K and (3^2 (Ladd and Howe |1997|) , a gas-to-dust ratio 
~ 100 and dust opacities of 0.005 cm^g^^ (Preibisch |1993| ), the masses associated with the 
two 222 GHz sources, within 1.4", are 2.5 Mq and 1.8 Mq. From the 110 GHz continuum. 



the mass associated with the envelope is found to be 13.6 Mq. Ladd & Howe (|1997D, for 



comparison, computed a mass of 10 Mq within a radius of 3000 AU. At 110 GHz, there is no 
spatial resolution to separate the two sources. 



7. IRAS 23011+6126 Spectral Energy Distribution 

The CVF ISOCAM observations clearly detect an embedded source that we identified 
with IRAS 23011+6126. Figure 12 shows the composite spectral energy distribution (SED) of 
IRAS 23011+6126. It resembles that of a Class I source (Lada |1987|) , rather than a Class 



as proposed by Lefioch et al ( |1996|) , based on millimetric observations. Presently, there are 
reasons to believe that the distinction between Class and Class I is not as sharp as previously 
thought (see e. g. Gregersen et al. pOOO|) . 



Following NCGM, model fits to the SED are calculated assuming: (1) a dust opacity 
dominated by bare silicates at a temperature of 18 K and a density of 7.5xl0^cm~'^, and (2) 
silicates with a thin ice mantle at a temperature of 18 K and a density of 6xl0^cm~^. The 
models assume power-law density and temperature distributions, with a core inner radius of 
0.065 AU and outer radius of 2500 AU. As Figure 12 shows, the continuum level of CAM 
observations near log(z/) ~ 13.3 is lower than that of the bare silicates model. As the success of 
the model using silicates with a thin ice mantle suggests, this can be attributed to absorption 
by H2O at 10-12 yum . The effect has been noticed in other low mass protostars such as Elias 29 



(Boogert |1999| ). This model produces a mass envelope of 13.2 Mq and bolometric luminosity 
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of ~ 34 Lq. Comparable masses have been obtained assuming a constant density over a 
spherical volume (10 Mq, Ladd & Howe |1997|) and in dust temperature dependent models 
(18 M0, Lefioch et al. |1996|) . The bolometric luminosity derived by our models, however, is 
a factor of 2 smaller that these cases ~ 30 L© versus 70 Lq or 100 L0. We suspect the 
discrepancy is due to the integration methods used, since coarse approximations overestimate 
the integrated flux. Our estimates of the luminosity and envelope mass are more appropriate 
for an intermediate mass object rather than a low mass Class I system; although not as 
massive and luminous as IRAS 05553+1631 (= GAL 192.16-03.82) (Shepherd et al. pJ98D or 



IRAS 20126+4104 (Shepherd et al. |2000D. 



8. Conclusions 

The Cep E outflow has been studied in the mid (~ 5 — 17 fim) and far infrared 
(~ 40 — 200 /im) using the ISO instruments CAM and LWS, and at millimetric wavelengths 
(110.2 and 222 GHz) using OVRO. In the near and mid-IR, the Cep E morphology is similar to 
that expected for a jet driven outflow, where the leading bow shocks entrain and accelerate the 



surrounding molecular gas (e. g. Raga & Cabrit |1993| ; Masson & Chernin |1993| ) . As expected, 
fine structure atomic/ionic emission lines are found, which very likely come from the leading 
bow shock, in both the Mach Disk and the stagnation tip, where J-shocks are dominant 
(Figure 13). The H2, H2O and CO molecular emission could arise further 'downstream' at the 
bow shock wings where the shocks are oblique and more likely to be C-type (Smith p!991| ). The 



range of shock velocities to excite these species is 8 — 35 km s~^. The ^^CO emission arises 
from entrained molecular gas and, as in other outflows, a compact high velocity emission is 
observed, together with an extended low velocity component that almost coincides spatially 
with the H2 near-IR emission. The millimetric continuum emission shows two sources. We 
identify one of them with IRAS 23011+6126 and we postulate that it is the driver of the 
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Cep E outflow; the other, also an embedded source, is likely to be driving one of other 
outflows observed in the region. Finally, the strong [C II] 158 /im emission must originate 
from an extended PDR, very likely excited by the nearby Cepheus 0B3 association. Our main 
conclusions are: 

1. The mid infrared spectra show strong pure rotational H2 lines, that coincide spatially 

with the H2 1-0 S(l) 2.12 fim emission. Excitation diagrams indicate that this emission 
is coming from a region at 1000—1300 K, in agreement with the post-shock temperatures 
expected for shock velocities in the range 20—30 km s^\ consistent with C-shock diagnostics 
of H2 line ratios. 

2. The far infrared spectrum is rich in H2O and CO lines. The line fluxes can be modeled 
using an LVG code, assuming that the H2O and CO emission arises from the same region. 
Two extreme models fit the observations, with temperatures either of 220 or 1200 K (Table 4). 
The resulting H2O abundances are between 100 and 4000 times greater than the one found 
in quiescent molecular clouds, confirming once more the presence of high water abundance 
in the shocked gas associated with young stellar outflows. From comparison with shock 
models, we conclude that C-shocks with Vshock ~ 8-14 km s~^ can reproduce the H2O and CO 
observations. 

3. Two strong atomic lines [O 1] 63 fim and [C 11] 158 fim are present in the LWS 

spectra. The ratio of their fluxes is close to 1, indicating a PDR origin. Taking into account 
the significant contribution to the [O I] line from J-shocks, and comparing with PDR models 
(Kaufman et al p.999| ), we find that two models reproduce the observations. Go = 30, n = 60 
cm~^ and Gq = 6, n = 5500 cm^^ erg s~^ cm~^, where Gq is the incident FUV flux in units of 
1.6 X 10~^ erg s~^ cm~^. Cep E is relatively close to the Cepheus 0B3 association, which can 
produce the diffuse FUV field required by these models. 

4. The bipolar outflow can be seen in the high resolution ^^CO image at 110.2 GHz. This 

emission is bounded by the H2 condensations. Outflow velocities extend from to -20 km s~^. 
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and are centered on the ambient cloud velocity of -11 km s~^. There is also evidence for high 
velocity gas at 30 and -40 km s~^ very close to the source. Prom the ^^CO observations we 
estimate, for the North and South outflows lobes respectively, time scales of 4000 and 8000 
years, masses of 0.08 and 0.05 M©, similar kinetic luminosities of O.IL© for both, and mass 
loss rates of 7 x 10~^ and 6 x 10"^ Lg/year. This confirms that the outflow is very young as 
the high water abundance seems to indicate. 

5. The continuum images at 222 GHz show the presence of two unresolved sources within 

the positional error bars of IRAS 23011-1-6126. The masses associated with them are 2.5 
Mq and 1.8 Mq. The lower spatial resolution of the continuum observations at 110 GHz allow 
us to detect the surrounding envelope which has a mass of 13.6 M© and a radius of 1100 
AU. No evidence is found for any of the other outflows observed at near-infrared of radio 
wavelengths. The presence of a second source, however, strongly suggests that the multiple 
outflows are due to multiple sources and not to precession. The separation of these sources is 
about 1.4"or ~ 1000 AU, and so even for an intermediate mass system, the binary period 
would be of the order of 7 x 10^ — 10^ yrs, which is longer than the kinematical age of the 
outflow. 

6. A simple envelope model for the spectral energy distribution of IRAS 23011-1-6126 

predicts a mass envelope of 13.2 and a luminosity of ~ 30 from silicates with thin ice 
mantles. The spectral features observed at 5— 17/Lim in IRAS 23011+6126 are closer to those 
of a Glass I source, and the mass and luminosity estimated are those of an intermediate mass 
YSO, rather than a low mass Glass object. 
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Table 1. OVRO observations 



Map 


Beam Size(") (PA°) 


RMS noise (l-a) 


110 GHz (cont.) 


3.1x2.8 (-86°) 


1 mJy/beam 


222 GHz (cont.) 


1.4x1.0 (-78°) 


7 mJy/beam 


C^^O (1-0) 


8x7 (+63°) 


0.05 Jy /beam/channel 


i^CO (1-0) 


8x7 (+54°) 


0.03 Jy /beam/channel 
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Table 2. Cep E LWS Integrated Fluxes^ 



A(/um) 


Element 


Transition 


North 


South 


63.18 


[0 I] 


3P1-3P2 


9.3(0.5) 


12.1(0.7) 


75.38 


0-H2O 


321-2l2 


1.8(0.4) 


0.9(0.2)*^ 


78.74 


0-H2O 


423-3l2 


0.6(0.5)^ 




79.36 


CO 


33-32 


1.3(0.4)"^ 




82.03 


0-H20 






0.9(0.4)^^ 


84.60 


OH 


3/2-3/2 (7/2-5/2+ 4-3) 


1.5(0.4)^ 


0.8(0.3)^^ 


89.99 


P-H2O 


322-2ll 




1.0(0.2)^ 


96.77 


CO 


27-26 


1.2(0.2) 


0.8(0.1) 


99.49 


0-H2O 


55-4l4 


1.2(0.4)^^ 


0.7(0.3)^ 


100.46 


CO 


26-25 


0.6(0.4)^ 


0.8(0.6)^ 


100.91 


0-H2O 


5 14-423 


1.0(0.4)<^ 




104.44 


CO 


25-24 


0.7(0.3) 




108.07 


0-H2O 


221-I10 


2.6(0.5)^ 


0.9(0.4)"^ 


108.76 


CO 


24-23 


1.0(0.8) 


1.1(0.6) 


113.54 


0-H20 


4l4-33 


2.7(0.3)f 


1.6(0.3)<^'f 


119.44 


OH 


3/2-3/2 (5/2+-3/2- 2-1) 


0.3(0.1) 


0.5(0.1) 


125.35 


P-H2O 


44-3l3 


0.3(0.1) 


0.4(0.2) 


126.71 


P-H2O 


331-322 


0.4(0.2) 


0.3(0.2) 


130.37 


CO 


20-19 


0.7(0.1) 


0.8(0. l)'^ 


137.20 


CO 


19-18 


0.7(0.1) 


1.0(0.2)-^ 


138.53 


P-H20 


3l3-22 


0.3(0.1) 


0.5(0.2)^ 


144.78 


CO 


18-17 


1.2(0.2)'^ 


1.4(0.4)^ 


145.52 


[0 I] 


3Po-3Pi 


<0.3'i 


0.9(0.6)^^ 


153.27 


CO 


17-16 


1.0(0.1) 


1.9(0.1)^ 
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Table 2 — Continued 





Element 


Transition 


North South 


157.74 


[C II] 


2P3/2-2P1/2 


7.0(0.3)'^ 8.1(0.2)^ 


162.81 


CO 


16-15 


1.7(0.1) 1.9(0.1)^ 


173.63 


CO 


15-14 


2.1(0.6)'^ 2.2(0.6)^ 


174.62 


0-H20 


33-2l2 


1.6(0.9)'^ 1.5(1. 0)<^ 


179.53 


0-H20 


2i2-li 


2.7(0.2)'^ 3.1(0.4)^ 


185.99 


CO 


14-13 


2.1(0.2) 2.3(0.2)^ 



'^In units of 10 ^^W/cm^. l-cr uncertainties are given 

in parenthesis. 

^FWHM was fixed at 0.29//m. 

^Blended with CO 29-28 90.16yum. 

<iFWHM was fixed at 0.60^m. 

^Blended with P-H2O 220-I11 100.98/xm. 

^Blended with CO 23-22 113.76//m. 

Horizontal dotted lines indicate non detections 



-30- 



Table 3. Rotational H2 0-0 Observed Fluxes'' 



Line(/xm) (15,12) (15,13) (15,14) (17,17) (17,18) (17,19) (19,21) 



S(7) 5.51 3.9(0.6 

S(6) 6.11 3.1(0.5 

S(5) 6.91 5.7(0.2 

S(4) 8.02 3.5(0.2 

S(3) 9.66 2.5(0.2 

S(2) 12.28 2.0(0.2 



4.3(0.5) 2.0(0.3) 

3.3(0.4) 1.0(0.2) 

7.2(0.3) 3.8(0.2) 

4.2(0.2) 2.0(0.2) 

4.3(0.2) 3.2(0.1) 

4.0(0.3) 3.1(0.2) 



2.3(0.3) 7.5(0.5) 

0.9(0.2) 4.1(0.2) 

5.1(0.3) 12.7(0.6) 

1.4(0.2) 2.6(0.2) 

1.9(0.1) 5.2(0.3) 

1.4(0.1) 1.9(0.2) 



4.1(0.2) 1.5(0.3) 

2.2(0.3) 1.5(0.3) 

6.5(0.3) 3.3(0.1) 

1.8(0.1) 1.6(0.1) 

3.8(0.2) 1.9(0.1) 

1.6(0.2) 1.3(0.2) 



'^Fluxes dereddened by E(B-V) = 1.0, in units of 10 W/cm^/sr 



Table 4. LVG parameters 



North South 



Temperature (K) 


220 


1200 


215 


1200 


n(H2) 


2.0 


0.02 


4.0 


0.04 


Wexp (km s"^) 


20-40 


25-40 


30-40 


15-60 


Source Diameter (") 


5.1-3.7 


5.7-4.3 


5.8-5.1 


9.5-4.7 


Nco (10^^-2) 


9.5-18.0 


16-26 


7.3-9.3 


2.8-11.3 


No-H20 (lO^W^) 


3-5 


30-50 


0.8-1.0 


5-20 


^p-H20 (lO^^cm-2) 


1-1.7 


10-16.7 


0.27-0.3 


1.7-6.7 


Leo (10-^^ W cm-2)^ 


2.1 


2.6 


2.4 


2.8 


K-H20 (10-^' W cm-2)b 


1.8 


2.1 


1.3 


2.1 


l^p-H20 (10-^^ W cm-2)- 


0.6 


0.6 


0.4 


0.6 



^Observed CO fluxes are 14.3(1.2) x 10"^^ (North) and 14.2(1.2) x 
10-19 W cm-2 (South) 

"^Observed 0-H2O fluxes are 14.2(1.4) x lO'^^ (North) and 
9.6(1.3) X 10-19 W cm-2 (South) 

^Observed P-H2O fluxes are 1(0.3) x lO-^^ (North) and 2.2(0.4) x 
10-19 W cm-2 (South) 
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Figure Captions 

Fig. l.~ ISOCAM-CVF observations of the Cep E outflow at H2 0-0 S(5) 6.1^m. The grey 
scale depicts the H2 0-0 S(5) 6.91 fj,m continuum-subtracted emission, and the superimposed 
contour image is the NIR vibrational H2 1-0 S(l) 2.121 //m hne (AYA). The FOV is ~ 2'. 
The locations of the LWS circular apertures are also shown in the image (r~ 40"), where the 
crosses represent the two pointing positions. 

Fig. 2. — ISO-LWS spectra of the Cep E North and South lobes after continuum subtraction 
and averaging over all detectors. Integrated line fluxes are given in Table 2. 

Fig. 3.— ISOCAM-CVF observations of Cep E outflow. In grey scale is the H2 0-0 S(5) 6.91/im 
image. Superimposed are the background subtracted spectra obtained at several positions along 
the outflow. Their scale is identical to that of the two magnified spectra shown in the upper 
left, for one position in the South and one in the North lobe. The spectrum at the central 
source position, IRAS 23011-1-6126, is shown in the lower right on a different scale. Below it, 
is the spectrum of the background averaged over 100 pixels. Orientation and pixel size are 
indicated in the figure. 

Fig. 4. — ISOCAM CVF images of the Cep E outflow at the rotational H2 0-0 lines, after the 
emission from the adjent line-free planes have been subtracted. All contour levels are from 0.3 
to 3.1 by 0.4, in units of 10~^ W/cm^/sr. The flgures have been plotted in different grey scales 
to enhance the contrast. The crosses indicate the position of the source. 

Fig. 5. — Comparison between the LVG predictions and the observations for CO and H2O lines 
fluxes, (a) CO line fluxes versus the rotational quantum number J„p for Cep E North lobe. 
Error bars indicate 1-a uncertainties. The two extreme LVG models that fits the observations 
are shown in solid (T = 220 K, n(H2) = 2x10^ cm'^ N(CO) = 8.65x10^^ cm'^ D = 5.1", 
Vexp = 20km s-i) and dashed (T = 1200 K, n(H2) = 2x10^ cm-^, N(CO) = 1.6x10^^ cm-^, D 
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= 5.7", Yexp — 25 km s ^) lines, (b) Same as (a) for Cep E South lobe. The two LVG models 
are in solid (T = 215 K, n(H2) = 4x10^ cm'^, N(CO) = 7.3x10^^ cm'^, D = 5.8", v^^p = 30 
km s"^) and dashed (T = 1200 K, n(H2) = 4x10^ cm'^, N(CO) = 1.1x10^^ cm'^, D = 15", y^xp 
— 50 km s~^) lines, (c) Ratio between LWS fluxes of 0-H2O hues and LVG model predictions 
for Cep E North lobe. Error bars represent 1-a uncertainties. Open squares correspond to 
the low temperature LVG model (same parameters as above with N(o-H20) = 3x10^'' cm^^), 
while crosses correspond to the high temperature model (with N(o-H20) — 3x10^® cm~^). (d) 
Same as (c) for Gep E South lobe. Open squares correspond to the low temperature LVG model 
(with N(o-H20) = 8x 10^^ cm~^), while crosses correspond to the high temperature model (with 
N(o-H20) = 1x10^^ cm-2) 

Fig. 6. — Determination of the H2 temperature and column density from the H2 0-0 rotational 
hues observed with ISOGAM-GVF (corrected by extinction, E(B-V)=1.0). The first panel is 
the dereddened observed spectra(sohd hue) at pixel (15,12), together with the model (dotted 
line) derived from the fitted temperature and column densities. On the right is the excitation 
diagram and the best-fit parameters for this pixel. The rest of the panels show the excitation 
diagrams for six positions along the flow with their corresponding best-fit values. 

Fig. 7. — GO/H2 (rotational) cooling ratio (grey scale) as a function of pre-shock density and 
shock velocity for G-shocks models (with B = (nj^/cm"^)^/^ jiG; from Kaufman and Neufeld 
1996). The coolings obtained from the 200 K and 1200 K LVG models for Gep E North and 
South lobes are shown in dashed and dotted lines respectively. Gontour levels are from to 6.4 
by 0.4. 

Fig. 8. — GO/H2O cooling ratio (grey scale) as a function of pre-shock density and shock 
velocity for G-shocks models (with B = (nij/cm~^)^/^ jiG; from Kaufman and Neufeld 1996). 
The coolings obtained from the 200 K and 1200 K LVG models for Gep E North and South 
lobes are shown in dashed and dotted lines respectively. Gontour levels are from -1 to 2.25 by 
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0.25. 

Fig. 9. — Superposition of the H2 1-0 S(l) 2.12 //m emission (in grey scale) on the ^^CO and 
222 GHz continuum in contours. The dashed contours towards the north correspond to the 
red-shifted lobe {vlsr = -9-5 to -2.5 km s~^), while the solid contours to the south correspond 
to the blue-shifted lobe ( vlsr — -20.5 to -13.5 km s~^). Contour levels for ^^CO are 0.4 to 3.0 
by 0.4 and 0.2 to 1.4 by 0.2 respectively (in Jy/beam). Contour levels for 222 GHz are from 
0.02 to 0.05 by 0.01 (in Jy/beam). The beam size is showed at the lower right corner. 

Fig. 10. — Position velocity diagram in the ^^CO J=l-0 transition along the axis of the Cep E 
outflow. The lowest contour correspond to a 2-a level. 

Fig. 11. — Left: same as Figure 9 but for the high velocity gas. The red lobe corresponds to 
the gas with velocities between and 20 km s~^ (dotted contours) and the blue lobe to gas 
with velocities between -47 and -20 km s~^ (solid contours). The contour levels are 10 and 15 in 
both cases. The detections are at 5-7 a level. Right: 222 GHz and 110 GHz continuum image 
centered on IRAS 23011+6126. The flux levels are from 0.02 to 0.05 by 0.01 and from 0.003 to 
0.02 by 0.003 respectively (in Jy/beam). 

Fig. 12. — Spectral energy distribution of IRAS 23011-1-6126. Data points are taken from Ladd 
& Howe 1997 (submillimeter; squares), from Noriega-Crespo et al. 1998 (IRAS and ISOCAM; 
open triangles and upper limit) and from our new OVRO results (filled symbols). The thick 
solid line represents the CVF spectra extracted at the source position. A model assuming a 
dust opacity dominated by bare silicates, at a temperature of 18 K and n = 7.5xl0^cm~^ is 
shown in sohd line. The dotted fine is for sihcates with a thin ice mantle at the 18 K and n = 
6xlOW-^ 

Fig. 13. — Schematic view of the different emitting regions in the Cep E outfiow adopting a 
jet-driven model. 
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Flux(W/cm //um/sr) vs. X(/xm) 
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Flux(W/cm /^m/sr) vs. A(^im) ln(N/g) vs. E (K) 
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